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I. INTRODUCTION 

Low dimensional integrable quantum systems, with or without open boundary conditions, which describe strongly 
correlated fermions M form an important class of lattice integrable models, which have attracted much international 
attention (see, e.g. Pk| and references therein). Recently, trying to extend the existing two component electron 
models to multi-component cases, we proposed Q an eight-state integrable model and its two-parameter (or q- 
deformed) version with Lie superalgebra gl(3\l) and quantum superalgebra E/ g [£/Z(3|l)] symmetries, respectively. One 
of the features of these two models is that they contain correlated single-particle and pair hoppings, uncorrelated 
triple-particle hopping and two- and three-particle on-site interactions. By eight- state, we mean that at a given 
lattice site j of the length L there are eight possible fermionic states: 



|0), ctjO), ct,|0), cUo), 



4i42l°>> ct lC t 3 |0> s cj. 2 ct 3 |0), ct lC t 2C t i3 |0), (LI) 



where cj a (cj, a ) denotes a fermionic creation (annihilation) operator which creates (annihilates) a fermion of species 

a = 1, 2, 3 at site j; these operators satisfy the anti-commutation relations given by {c\ a , Cj^p} — Sij5 a p. 

Recently we formulated in j?J a general and fully supersymmetric version of the boundary inverse scattering method 
HP) , and constructed a large number of integrable boundary conditions |T(| for various models of strongly correlated 
fermions. In this paper, we continue our study of open boundary conditions and consider the integrable eight-state 
fermion model with C/ 9 [<7/(3|l)] symmetry. We present nine classes of boundary conditions for this model, all of which 
are shown to be integrable by the graded boundary QISM (?J . We solve the boundary systems by using the coordinate 
Bethe ansatz method and derive the Bethe ansatz equations for all nine cases. 

II. OPEN BOUNDARY CONDITIONS 

We consider the following Hamiltonian with boundary terms 

L-l 

H = H jd+i (g, «) + H* oundaiy + ff r b ° undary , (II. 1 ) 

3=1 
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where ^/^ oundar y ; jj ° nn aT y are respectively left and right boundary terms whose explicit forms are given below, and 
Hjj+i is the Hamiltonian density of the two-parameter eight-state supersymmetric fermion model || 

Hj,j+i(g,K) = -^(ct Q Cj+i, Q +h.c.)exp I ~(t} + k) ^ n j+e(p _ a)i0 - -(f) - k) ^ n j+1 _ e{fJ _ a)Jj 

C / \ I sinhft ^ / t t in 

+ 2 2^ \ n hP n 3,i + %-+i,^%+i,7J f ~ 2sinhKfo + 1) ^ ( c ].» c 3 ,p c j+iJJ c 3+i,a +h.c.) 

cxp |-(| - sign( 7 - 2)/c)n i>7 - (~ + sign(7 - 2)«)n i+ i jT | 

2 cosh k sinh 2 k / + + + \ 

' sinh K (. 9 + l)sinh K ( g + 2) (^2^+1,3^+1,2^1,1 + h.c. ) 

tea -ko sinh/i v 

+e 5 nj + e 9 n j+ i - 2s i n Y 1K ( g + i) 2^ W^n^ + n j+1>a n j+1 ,p) 

2cosh/t(g + 1) sinh 2 k , , ,_ „ s 

+ sinh.( g +l)sinh.( g + 2) ( ^ in ^- 3 + ' 1 • ' ' ' ' ( IL2 ) 

where g, k are two free parameters, nj = + rij t 2 + %,3 with rij ta = c| Q Cj jQ being the number operator for the 
electron of species a at site j, 6(f3 — a) is a step function of (/? — a) and 

, = _ h *^L_, C= lln . ^ 2 "(g+l) C =-ln . f^^ * (n-3) 

sinhft(<7 + l) 2 sinh ng sinh K(g + 2) sinh n(g + 2) 

As is shown in ||, the symmetry algebra underlying (EL 2) is quantum superalgebra t/q[ g Z(3|l)]. The parameter k 
is related to the deformed parameter q by q — e K . 



i related to tne deiormed parameter q by q = e . 
We propose the following nine classes of boundary conditions: 

n ,.s ^boundary sillh K g / _ sinh K 

Case (i) : if K - y = ^— e 5 ni _ — (ni^rai^ + ni^n^a + rii,ini, 3 ) 

sinh -2~K V sinh(l + 



2 f£ 
2 sinh kcos!i(1 + ^-)k 

H 77 77 ™1,1™1,2™1,3 I , 

sinh(l + -y)k sinh(2 + -^-)k 



^boundary sinh Kg / f+ sinh K 

sinh^-K y sinh(l + -£-)k 

2 sinh 2 KCOsh(l + %)k \ 

H 77 77 — «l,i«l,2^l,3 ; (II.4) 

sinh(l + ^-)Ksinh(2+ y 

n / • • \ ^boundary sillh K g / _l" K , . Smh K 

Case (11) : if H y = ^— e 2 (ni,2 +711,3) 777 — "1,2x11,3 

sinh-^-K y sinh(l H — i-)«; 

^boundary sillh K ff | f Smh K \ 

H n = 777- e 2 (ni,,2 + ri L , 3 ) 7J7 — n Lt2 n Lt3 \ ; (II.5) 

sinh-^-K y sinh(l + -^-)k 

Case (iii) : ^oundary = sinh/tg e _£^i K 



sinh k 



/// 

"1,3, 



rrboundaty _ sinh Kg -t— K 



= ~ 7777-e 2 K n L x, (11.6) 

sinh — k 
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Case (iv) : ^boundary = I e -^ m — (ni,im, 2 + ni, 2 ni, 3 + m.im.a) 

sinh — k y sinh(l + — )k 

2 sinh 2 k cosh( 1 + %■ ) /c 

H 77 77 "1,1711,2X11,3 

sinh(l + -2-)«;sinh(2 + -^-)k 
^boundary sinh Kg / «|* sinh k t 

#rt = — e ^ K (n L , 2 + "L, 3 ) 777 — n L . 2 n L . 3 | ; (II.7) 

sinh-5-K V sinh(l + -^-)k 



sinh Kg / _£^„, sinh/t 



Case (v) : ^boundary = ~g e -^ K(m 2 + ni a) _ 



sinh-^-/t y sinh(l + -^-)k 



ni,2ni,3 



^boundary sinh Kg I f± Smh K 

# rt = 77— e 2 «l 77 — (™l,i"l,2 + n La n L ^ + n L ,in L ,3) 



sinh^K y sinh(l + -£-)k 



2 £ 7 
2sinh kcosIi(1 + ^-)k . 

H 77 77 n L,l n L,2 n L.3 (H. 

sinh(l + ^-)Ksinh(2+ / 



n / -\ ^boundary sinh Kg / _£- K sinh K 

Case (vi) : H lt = — e ^ m -j — (rai,im,2 + ni, 2 ni, 3 + ni,mi, 3 ) 

sinh-5-K y sinh(l + - 2 z -)k 

2 sinh 2 k cosh( 1 + % )/t 



t/f n i, l n l, 2 n l, 3 

sinh(l + -j l )ac sinh(2 + -^-)k 

rrboundary Smh Kg £+_ 

sinh — ^ — k 



Case (vii) : ^boundary = sinh k. 9 e _!^L K 



£iu - "1,3, 

sinh — k 



^boundary sinh «0 / £± K sinh K , , 

ff rt = 77— e 2 n L 77 — (ni,in i: 2 + n L ,2n L ,3 + n L ,\n L fi) 

sinh-^-K y sinh(l + -4-)k 



2 sinh 2 k cosh(l + %)k 



H 77 ^7? — n LA n Lj2 n L>3 ; 

sinh(l + ^-)«sinh(2 + y 

Case (viii) : = f e"->i,2 + m, 3 ) ^77-^^1,3 ' 

sinh—K y sinh(l + —)k 



(11.10) 



tjboundary Smh K<7 1±_ K 

= 7777— e 2 (II.ll) 

sinh — v — k 

Case (ix) : ^oundary = sinh«g e -l^„ nii3i 

sinh — — k 

^boundary sinh/«g [ f±. sinh k \ mio > 

= 777— e 2 (n L ,2 +«L,3j 777 «L,2«L,3 5 (11.12) 

sinh-^-K y sinh(l + ^-)k y 

where £,±(a = I, II, III) are some parameters describing the boundary effects. As will be shown in next section, all 
nine classes of boundary conditions are integrable. 
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III. BOUNDARY K-MATRICES AND QUANTUM INTEGRABILITY 



Quantum integrability of the boundary conditions ( |n.4| - [lLT2 ) can be established by means of the (graded) boundary 
QISM recently formulated in m . We first search for boundary K-matrices which satisfy the graded reflection equations: 



1 2 2 1 

R\2{u\ - u 2 ) K- (ui)R 2 i(ui + u 2 ) K- (u 2 ) =K- (u 2 )Ri 2 {ui + u 2 ) K- (ui)i?2i(wi - "2), 
R s £ tst2 (- Ul + u 2 ) K* 1 {u x )R 12 {-u x -u 2 +4) Kf> (« a ) 

=K l f 2 (u 2 )R 21 (- Ul -u 2 + A) (u^R^i-U! + u 2 ), 



(III.1) 



(m.2) 



where R(u) £ End(V <g> V), with V being 8-dimensional representation of C/ 9 [t//(3|l)], is the R-matrix of the two- 
parameter eight-state supersymmetric fermion model |5), and R 2 \(u) — P\ 2 Ri 2 (u)P\ 2 with P being the graded 
permutation operator; the supertransposition st^ (/x = 1, 2) is only carried out in the /i-th factor superspace of V ® V , 
whereas ist^ denotes the inverse operation of st^. 

The whole procedure of solving the reflection equations is quite involved. We shall not spell out the details, but 
state that there are three different diagonal boundary K-matrices, KL{u), KL^u), X£ /7 (u), which solve the first 
reflection equation ([ILl): 



Kl(u) 



Kt! (u) 



K ul {u) 



where 



£L ££. ii 

sinh — K sinh(l + -f-)K sinh(2 + -=-)n 



j — diag (Ai (u) , BL (u) , BL (u) , B J _ (u) , Ci (u) , C£ (u) , C£(m), D[_ (u)) 



1 



sinh -=-fcsinh(l + )i; 
1 



—diag (A£'(u), AUiu), B'Jiu), B£'( U ), B'J (u), BL 1 (u),C H (u),C H (u)) 



sinh — — k 



— diag ( A 111 («) , A 1 ! 1 (u) , A 111 (u) , B ni (u) , A 111 («) , B ni («) , B ni («) , b£" (u)) , 



A£(u 
B£(u 
C£(u 
L>£(u 
A£ 7 ( u 

c£ 7 ( M 

A^iu 



-e 2 sinh Ksmh — 



u . — 4 — ££ -|- u 
— Ksinh k, 



inn • v. £ + « . , -2 - ££ + « . , -4 - ££ 

e 2 sinh — - — k sinh k sinh — 



_i,. K . ££+M . 2 + ££+it . -4-££+u 

— e 2 sinh Ksinh Ksinh k. 

2 2 2 

_3 MK .,££+" 2 + C£+u . ,4 + f£+u 
e 2 sinn k sinn k smn k, 



e UK sinh 
— sinh 



2 2 
-££ J + u . -2 - ££ J 



-k sinh ■ 



2 2 
££ J + u . -2 - C£ 7 + u 



sinh 



= -e 2K sinh- 



2 

e£ J + u 
2 



-K sinh ■ 



k sinn k, 



111 



-K, 



B HI {u) = e~* K sinh 



(III.3) 



(III.4) 



The corresponding K-matrices, if£(u), K¥(u), K+ (u), can be obtained from the isomorphism of the two reflection 
equations. Indeed, given a solution Kt(u) of (III.l), then Ki_{u) defined by 



Kl s \u) = MK°l(-u + 2), a = I, II, III, 



(111.5) 



are solutions of (III. 2). The proof follows from some algebraic computations upon substituting ( III.5| ) into ( III.2| ) and 
making use of the properties of the R-matrix . Here M is the so-called crossing matrix, which is given by in the 
present case, 



M = diag (1, 1, e 2K , e 4K , e 2K , e 4K , e 6K , e 6h 



(HI.6) 
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Therefore, one may choose the boundary matrices K^_{u) as 

Ki(u) = diag (4+.(u), B' + {u), e 2fi Bl(u), e 4K .B|(u), (u), e 2K C|(u), e 4K C|(u), , 
K z + Z (u) = diag (^( U ), A%(u), Bi^u), ^(u), B^(u),C"(u),D^(u),D^(u)) , 



Kl u (u) = diag (A I + II (u),A I + II (u), e a, M£"(u), e 2,t ^"(u), B£"(u), e 2K B+ // (u), e 2K B" / (u)) , (III.7) 



where 



CI 



.4 



S|' 



ci 1 



5 



+ 



2.9 - £ 



«) 



7i) = e -(i" +2 ) B sinh 2g ^+ 

u)=e (i«-2)K sinh ^_i4_l 

; 2 



u 2 ff - 2 - £| 
— k sinn — 



— k smn — 



u 2g - 4 - e+ 
— K smn 



- k sinh 



2 



u 2g + 2 - £| 
— k smn 



— k smn — 



-k, 



u) 
u) 



e 2 " K sinh 



2.g - e( 



sinh 



2 

2.9 -& 



u 2g + 2 - C| 
— K smn 



u 2 5 + 4 - C| 
— K smn 



-K, 



u . 2 g - 2 - e| 7 

— k smn — ! — 



m) = e K sinh 



2.9 - 2 - e 



u) = e 2K sinh 



2 



u 2 ff + 2 - e| 7 

— K smn — 



-K, 

■</. 



2 

■« . 2, 9 + 2 - e| J 

— re smn 



-K, 



i) = e {u+2)K sinh 



2.9 + 2 - £ 



1/ 



u 2 S + 4 - e 

— k smn 



+ 



u) 
u) 



sinh 



2.9 -e 



III 



-K, 



e(^ S inh 2g + 4 -^ 



(III.8) 
(III.9) 

(111.10) 
(III.ll) 

then it can be shown that [t(ui), t(-u 2 )] = 0. Since K±{u) can be taken as K±(u), K±{u) and K± n (u), respectively, 
we have nine possible choices of the boundary transfer matrices: 

T (a ' b \u) = str (Kl(u)T(u)K h _(u)T- 1 (-u)) , a,b=I, II, III, (111.12) 

which reflects the fact that the boundary conditions on the left end and on the right end of the open lattice chain are 
independent. 

Now it can be shown that Hamiltonians corresponding to all nine boundary conditions are related to the second 
derivative of the boundary transfer matrix T^ a,b \u) (up to an unimportant additive constant) 

2sinhK g H(a b) 



Following Sklyanin's arguments ||] , one may show that the quantity 71 (u) given by 
T_(«) = T{u)K_(u)T- 1 {-u), T(u) = R 0L {u)---R 01 {u), 
satisfies the same relation as K_{u): 

1 2 2 1 

R 12 (ux - u 2 ) T- {u 1 )R 21 (u 1 + u 2 ) T- {u 2 ) =T- {u 2 )Ri 2 {ui + u 2 ) T- {ui)R 2 i(u 1 - u 2 ). 
Thus if one defines the boundary transfer matrix t(u) as 

t(u) = str(K+(u)T^(u)) = str (K+(u)T(u)K-(u)T- 1 (-u)) , 



H 

jj(a,b) 



K 

_(«,&)" 



(0) 



L-l 



4(V + 2W) 



i 



2(V + 2W) 



str K a + (O)Glq 



+2 str K a '+ {0)H? + str [K a + (0) (H 



l L0) 



(III. 13) 
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where 



Hi 



V = strvK% (0), W — str a 
+1 ~ Gj.j^ 



o 

K a . 



(0)H 



I',., ■ : , . i ((•)• 



(III. 14) 



Here Pj,j+i denotes the graded permutation operator acting on the j'-th and j + 1-th quantum spaces. ([II. 13) implies 



that the boundary two-parameter eight-state supersymmetric models admit an infinite number of conservation currents 
which are in involution with each other, thus assuring their quantum integrability. 



IV. COORDINATE BETHE ANSATZ ANALYSIS 



Having established the quantum integrability of the boundary models, we now solve them by using the coordinate 
space Bethe ansatz method. Following p,|l0|,f7[ , we assume that the eigenfunction of Hamiltonian (II. 2) takes the form 



1*)= iSf o ll ,-,a N {xi,---,X N )cl iai ---cl NaN \0), 

{(Zj.Oj)} 

N 

^ au --,a N (xi,- ■■,x N )=Y2 e P A aQu .... aQN (k PQ1 , • ■ • , k PQN ) exp(i ^ k Pj Xj), (IV. 1) 

p 3 =i 

where the summation is taken over all permutations and negations of k\ , ■ ■ ■ , k^ , and Q is the permutation of the N 
particles such that 1 < xqi < ■ ■ ■ < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the eigenvalue equation H\^>) = E\^S), one gets 

-A. ■ ■ .ctj .at,- ■ ■ (' ' ' : kj : kii ' ' ') — Sij {kii kj)A... aia ^...( K - • • , ki , kj ; , • ■ •) , 

Aon,--- {— kj ,■■■) — s L (kj; pi ai )A ai ...(kj, ■ ■ •), 

A..., ai {- ■ • , -kj) = s R (kj;p Lai )A... tai (- ■ ■ , kj), (IV.2) 
where Sij(ki,kj) are the two-particle scattering matrices, 
Sij(ki, kj)®^ = 1, a = 1,2, 3, 

Stf(fci,fcj)g6= • 5 ™ {X \ Xj \ a^b, a, b =1,2,3, 
sm(A, — Xj — in) 

Sij (fci, kj)t = e ^°(»-")^-^) Sm [ K 0,6=1,2,3, (IV.3) 

sin(Aj — Xj — in) 

where Xj are suitable particle rapidities related to the quasi-momenta kj of the electrons by 

fc(A) = 2 arctan(coth c tan A) , (IV.4) 

where the parameter c is defined by 

If sinh|(rj + «;) 1 (IV 5) 

4 \ sinh h(v ~ K ) i 

s L (kj\p\ ai ) and s R (kj\pL ai ) are the boundary scattering matrices, 

1 -piQ.e 1 ^ 



s (kj-,pi ai ) 



1 — PlQi<3 lkj 
1 I'l.n ( ' 



s R (k j;PLai ) = i PLai \ k . e*M*+D (IV.6) 
with p\ ai and pL ai being given by the following formulae, corresponding to the nine cases, respectively, 
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. _ KB __L K sinh Kg 

Case 1 : pi,i = pi,2 = Pi,3 = Pi = -e 9 + e 2 



sinh — k 



Pl,i = PL, 2 = PL, 3 =Pl = -e K9 + e 



f± K sinh Kg 



sinh -^k 



Case ii : pi,i = — e Kg , Pi,2 = Pi,3 = — e K9 + e 



PL,i = ~e K9 , PL,2=PL,3 = -e K9 + e 



e~ Kg + e" 



f- „ sinh Kg 
•■±- K sinh Kg 



Case hi : p M = pi j2 = -e~ K9 , pi i3 

/>/..: /'/...! - ' 7i . /'/.a • 
Case iv : p M = p 1>2 = Pi, 3 = Pi 



sinh -5-K 
sinh Kg 



sinh — 2 — ^ 
{ + K sinliKo 

K 9 _|_ c — 5 — K J 

sinh — ^ — k 



-ko , - — « sinh Kg 

e 9 + e 2 



sinh — k 

/>/..: ' ' /'/..J /'/..:> ' ' ''' • ' ~" 



Casev: p 11 = -e Kg , p 12 =pi_ = -e Kg + e 



C 11 ' 
sinh -5-K 

_ K sinh Kg 

sinh — k 



_„. f±„ sinhKo 
PL.l = PL,2 = PL,3 =p L = -e Kg + e 2 

sinh -^k 

. _ Kn _lL K sinh Kg 
Case vi : pi,i = pi, 2 = Pi,3 = Pi = -e g + e 2 ) 



sinh — k 



- ,. <•■</ . .— ■ - illliK // . 

^1 11 i 

sinh — ^ — k 



PL,i=Pi,2 = -e K9 , Pl,3 = -e K9 + e 2 



€ — Slllll 

Case vii : p ltl = p 1>2 — -e~ Kg , p 1 ^ = -e~ Kg + e ~ K ' ' 



sinh — 2— k 



Pl,i = PL, 2 = PL, 3 = PL = -e K9 + e 



f+ K sinh Kg 



sinh -£-K 



Case viii : pi,i = — e Kg , pi,2 = Pi,3 = — e Kg + e 



PL.\ = PL, 2 = -e a , Pl, 3 = -e 



Case ix : pi,i=pi, 2 = -e K9 , p 13 = -e K9 + e 



_ K sinh Kg 
sinh — ^— k 
sinh Kg 



sinh — ^ — k 



{ " J sinh Kg 



sinh — 2 — ^ 



PL,i = -e K9 , _pl,2 = PL,3 = -e K9 



s + „ sinh kg 



sinh -±-k 



(IV.7) 



(IV.8) 



(IV.9) 



(IV.10) 



(IV.ll) 



(IV.12) 



(IV.13) 



(IV. 14) 



(IV.15) 



As is seen above, the two-particle S-matrix (IV. 3) is nothing but the R- matrix ol the ?7,j[gl(3)]-invariant Hciscnberg 
magnetic chain and thus satisfies the quantum Yang-Baxter equation (QYBE), 



Sijiki, kj)Su(ki, ki)Sji(kj, ki) — Sji(kj, ki)Su(ki, ki)Sij(ki, kj). 



(IV.16) 
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8 



s R (k J ; PLai )^e^ L+1 h(k J ;p LA ) 

/l 

V 



Case ix : s L (kj\p\ ai ) = s{—kj]p\^\) 




(IV.26) 



1 

o 



s R (k 3 ;p Lat ) = e^ L+1 ^(k 3 ; PL . 1 ) 



( 1 

e 2 ^ S ? n ^ + T^ 

sin(C + +Aj) 






_2iA, sin(C + -Aj) 



(IV.27) 



V C ' Bi„(C + + A. ; ) 

Here / stands for 3x3 identity matrix and pi+, pl+ are the ones given in ( [V.8D ;(±, C± are parameters defined by 



9 ?± >/ 9 ?± 



2/ 



2/ 



(IV.28) 



We immediately see that (IV. f9) are the trivial solutions of the reflection equations (IV. 17), whereas (IV. 20) and 
(IV. 21) are the diagonal solutions |^,||. 

The diagonalization of Hamiltonian (II. 2) reduces to solving the following matrix eigenvalue equation 

Tjt = t, j = l,---,N, (IV.29) 

where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

Tj = Sr(fc^(_^;p lffJ )iij(fc i )i^(fc J a fl (fc J sp i a J )£5 h (fci) (IV.30) 

with 



Sf(kj) = S jtN {kj,k N ) ■ ■ ■ S jtj+ i(k j ,k j+1 ), 
Sj (kj) = Sj,j-i{kj, kj-i) ■ ■ ■ Sj^i(kj, ki), 
R j(kj) = Sij(k u -kj) ■ ■ ■ Sj-uikj-u-kj), 
Rf(kj) = Sj + i,j(k j+ i,-kj) ■ ■ ■ S Nd (k N , -kj). 



(IV.31) 



This problem can be solved using the algebraic Bethe ansatz method. The Bethe ansatz equations for all the nine 
cases are 



Mi 



e i^(L + i )F{k .. pi+:PL+) = yi 



sin(Aj - A^ + in/2) sin(Xj + A^' + in/2) 
sin(Aj - A^ - in/2) sm(Xj + A^ - in/2) ' 
G(A(1) . c c+) Q sin(A^ - A< 1] + in) sin(Ai 1} - A a) + in) 



(i) 



V sin(Ai 1} - Xj + in/2) sin(A^ 1} + Xj + in/2) 



n 



i sin(Ai 1} - Xj - in/2) sir^A^ + Xj - in/2) 

M 2 

n 

P =i 



sin(A a) - A y p> - in/2) sin(A^ + A)?' - in/2) 



(2) 



(1) A 00 

p 

'1) J_ A'---' 



p = i sin(A^. > - Ay - in) sin(A. 



0)_ A 0) 
tip 



in) 

a = !,•••, Mi, 



Mi 

n 

P =i 



j sin(Ay ) - A p iJ + i«/2) sin(A^ J + A p 1} + ire/2) 
sin(A 



(2) 



A 



(i) 



in/2) sin(AL 2) 7 + Lp 1 ' - in/2) 



(i) 



sin(AS 1) - A [ p ] + in/2) sin(Ai 1) + A ( p 2) + in/2) 

K(KW-t> c> \ TT sin(A ^ - A ^ ± ^ sin(A ^ ± A ^ + ^ 
1 ' U ' ^ ii sin(A< 2) - A p 2) - i«) sin(A^ + A p 2) - i«) ' 



P#7 



7 = M a , 



(IV.32) 



where 



F(kj-,p ll i,pL,i) = s(kj;px t i)s(kj;pL,i), (for all cases) , 
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G(AW;C-,C4 



sin(C-+A< r 1) + ^) S in(C++A^> + ^) 
sin(C--A< r 1) + ^) sin(C + -A< r 1) +^) 



S in(C + +A^' + ^) c2iA q) 

sin(C + -A< r 1) + ^) 

sin(C-+A< r 1 ' + ^) 2jA q) 

S in(C_-A< r 1) + ^) 

1 

1 

sin(C_ + A( r 1 ' + ^) e _ 2iA a) 
sin(C--A< r 1) + ^) 
S in(C + +A^' + ^) c2iA q) 
sin(C + -Ai 1) + ^) 

1 
1 

sin(C^+A< f 2) +z K ) sin(C++Ay'+ZK) 
sm(C^-A<, 2) +iK) sin(C^-A<, 2) +iK) 



n((;+Af+«) 2*A<, 2 > 



n(C;-A<, 2) +i K )^ 
n(C-+Ay ) +»«) c -2iA( 2 ^ 



-A<, 2) +i K )' 



n(CL 

n(C;+A^ 2) +m) c2iA (2) 



n(C;-A<, 2) +i K ) 
n(Cl+A<, 2) +»K) c „ 2 ^ A (2) 



n(C -Ai, - 1 +ik) 



case (i) 
case (u) 
case (in) 
case (w) 

case (u) 

case (w) 
case (vii) 

case (viii) 

case (ix) 

case (i) 
case (ii) 

case (Hi) 

case (iv) 
case (v) 

case (vi) 
case (vii) 
case (viii) 
case (ix) 



(IV.33) 



The energy eigenvalue E of the model is given by E 
coming from the chemical potential term). 



2 X^fLi cos fej (modular an unimportant additive constant 
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